Generalization of the 
Berezinskii-Kosterlitz-Thouless theory to higher 

vortex densities 

Carsten Timm 
Universitat Hamburg, 1. Institut fiir Theoretische Physik, 
Jungiusstrafie 9, D-20355 Hamburg, Germany 

February 23, 1996 

Abstract 

The Berezinskii-Kosterlitz-Thouless theory for superfluid films is gen- 
eralized in a straightforward way that (a) corrects for overlapping vortex- 
antivortex pairs at high pair density and (b) utilizes a dielectric approx- 
imation for the polarization of the vortex system and a local field cor- 
rection. Generalized Kosterlitz equations are derived, containing higher 
order terms, which are compared with earlier predictions. These terms 
cause the total pair density to remain finite for temperatures above the 
transition so that it is not necessary to introduce an ad hoc cut-off, as 
opposed to the original Berezinskii-Kosterlitz-Thouless theory. The low- 
temperature bound pair phase is destabilized for small vortex core energy. 
The behaviour of the stiffness constant and of the correlation length close 
to the transition is not affected by the higher order terms. A first-order 
transition as suggested by other authors is not found for any values of the 
parameters. The pair density is calculated for temperatures below and 
above the transition. Possible experiments and the applicability of the 
extended approach are discussed. The approach is found to be applicable 
even in a significant temperature range above the transition. 



1 Introduction 

Two-dimensional and quasi-two-dimensional systems with a two-component or- 
der parameter and short-range interactions typically fall into the universality 
class of the two-dimensional XY model. The behaviour of such systems, e.g., 
superconducting films, superfluid films and thin planar magnets, close to their 
phase transition is usually well described by the Berezinskii-Kosterlitz-Thouless 
(BKT) theory |. The discovery of the high-temperature superconductors. 
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which consist of weakly coupled superconducting layers, has led to renewed in- 
terest in this theory |3| . In BKT theory the phase transition is supposed to take 
place through the unbinding of spontaneously created vortex-antivortex pairs. 
For temperatures up to the transition temperature the predictions of BKT 
theory typically agree quite well with experimental results. There are several es- 
sentially equivalent mathematical formulations of the BKT theory ^, ||, 0| . 
Here, we follow Halperin's presentation j^, which is founded upon the heuristic 
approach of Kosterlitz and Thouless |^ . 

The BKT theory treats the interaction between different vortex-antivortex 
pairs through a suitable renormalization of the intra-pair interaction. This 
approximation is best justified if pairs are small and far apart. The so-called 
Kosterlitz equations, which are derived in BKT theory, relate the density of 
vortex-antivortex pairs of given size and the screening of the vortex-antivortex 
interaction at given separation to one another. Nelson and Kosterlitz and 
Halperin predict the general form of higher order terms of these equations, 
but do not derive the coefficients. Amit et al. utilize field-theoretical 

methods to derive a universal relation between the coefficients of two higher 
order terms. One would like to know the higher order terms, including the 
coefficients, to see if the typical BKT behaviour, which is found in experiment, 
survives in their presence. 

Without these terms the BKT theory suffers from the following weakness: 
If taken literally it predicts a diverging vortex density for T > Tc- This result 
cannot be correct. The creation of one vortex-antivortex pair costs at least 
twice the core energy. The density of pairs should thus be governed by a Boltz- 
mann factor containing this energy. In the BKT theory, a cut-off length ^_|_ is 
introduced ad hoc to obtain finite results. The length ^+ is meant to represent 
the average separation between free vortices. In a system of zero overall vor- 
ticity, all vortices and antivortices can be grouped into pairs. Thus f+ is not 
very well-defined. The concept of "free vortices" does not naturally fit into the 
theory. It is unsatisfactory that such a cut-off is necessary at all. It seems that 
an important contribution to the density has been neglected. 

Different objections arc put forward by Gabay and Kapitulnik [11|. These 
authors argue that one should use the Clausius-Mossotti formula [12| for the 
dielectric constant to take local field effects into account. Several authors 
[|3[ |l^ |l^ predict the breakdown of the BKT theory for small vortex core 
energy and the appearance of a first-order transition. These papers are inspired 
by Monte Carlo simulations |l^ , which indicate the presence of such a tran- 
sition for low vortex core energies. Minnhagen and Wallin employ a linear 
screening theory to obtain generalized Kosterlitz equations. It seems that local 
field effects are not included in their approach (see below) and that the restric- 
tion that vortices may only be created in neutral pairs is not taken into account. 
Furthermore, the authors use different sets of self-consistency equations below 
and above the critical temperature, whereas one would welcome an approach 
which predicts a phase transition without having to assume its presence be- 
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forehand. In Refs. [Q and ||T^ phase diagrams are derived that do not match 
the BKT resuh for high core energies, Ec > q^, where the BKT theory should 
be apphcable. Here, is the couphng constant appearing in the unscreened 
vortex-antivortex interaction V(r) = (7^1n(r/ro). Experimental evidence for a 
first-order transition does not seem to exist as of yet. 

In the following sections, a generalized BKT theory is proposed, which is 
based on straightforward physical ideas and takes the points mentioned above 
into consideration. A geometric correction, which avoids the diverging density, 
is employed. For the derivation of the screening of the intra-pair interaction, 
we start from an exact expression of linear response theory and then take the 
polarization of the vortex system and local field effects into account. We end 
up with higher order terms in both Kosterlitz equations. These terms are com- 
pared with the predictions mentioned above and the generalized equations are 
discussed. It is shown that these equations describe the vortex system even in 
a significant temperature range above the transition. 



2 Geometric correction at high pair density 

Let t'(r) dr be the areal density of vortex-antivortex pairs of sizes between r and 
r -I- dr. Then the total density of pairs is 

dri>{r), (1) 

where ro is the minimum pair size, which is of the order of the Ginzburg-Landau 
coherence length. 

The BKT ansatz for the pair density is 

;y(r) dr^^ 2TTr dr exp(-2/3Se) exp[-(3Ves{r)], (2) 

^0 

where A^o is the number of possibilities to place a vortex in an area r^, f3 = 1/T 
is the inverse temperature (we set fc^ = 1 throughout the paper), Ec is the core 
energy of a vortex and Vcs is the effective intra-pair interaction. The expression 
for the pair density contains a geometric factor 27rr dr, which represents the 
area of the circular ring in which the separation vector resides, and a Boltzmann 
factor exp(— /3[2£'c -I- K!ff(r)]), which contains the energy of the pair. 

We now understand why the total density diverges for T > Tc- In this tem- 
perature regime we expect the effective interaction to be screened, i.e. limr^oo Kiff (j") = 
const. Thus v{r) dr oc rdr for large r. Equation then implies that the total 
density n diverges with the system size squared. As noted above, this result can- 
not be correct since the energy of a pair is at least 2Ec > 0. Hence, an important 
effect must have been neglected in Eq. (||). As explained in the following, the 
error lies in the enumeration of pairs. 
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The idea of BKT theory is to group the vortices and antivortices into neu- 
tral pairs and to describe the interaction between different pairs through the 
renormahzation of the intra-pair interaction. However, we must not group the 
vortices and antivortices arbitrarily. For the approximation to be best justified, 
we must form pairs as small as possible We now imagine that we fill the 
system with pairs, starting from the smallest ones (of size tq). The probability 
that we add a pair of size r contains the geometric term and the Boltzmann 
factor present in Eq. (H). Furthermore, we must make sure that the enumera- 
tion remains consistent. If we placed the vortex of the new pair closer than a 
distance r' to the antivortex of another pair, where r' is the size of the latter 
pair, the new vortex and the old antivortex should be considered a pair. In 
that case, we would have miscounted the pairs. Our ansatz is to introduce an 
additional factor in the density, which denotes the probability that the pairs 
are placed properly. This probability is determined by the distribution of the 
smaller pairs already present. 

A new pair of size r be tentatively placed anywhere in the system. The 
probability that the new vortex is not placed closer than r' to the antivortex of 
any pair of given size r' is 1 — vir') dr' 7rr'^, since J^(r') dr' is the density of such 
pairs. The probability that the new vortex is not placed closer than r' to the 
antivortex of any pair of size r' < r is then given by 

[| [1 - i^ir') dr' 7rr'2] = exp [-^{r') dr' nr'^] 

(3) 



r' <r 



exp 



The probability that the placement of the new pair does not violate the enu- 
meration rule is just the square of the last expression since a corresponding 
condition must be satisfied for the antivortex. Thus we obtain 

,y(r) dr = 27rr dr exp{-2/3Ec) exp\-- pV^sir)] 

^0 



X exp 



-277 / dr'r'2jy(r') 



(4) 



instead of Eq. (||). 

This integral equation for ^{r) can be solved by transformation into a dif- 
ferential equation for Y{r) — J^^ dr' r'^v{r'). The solution is 

^ 2'Krn* exp [-(3V^s{r)] 

l + ^Ti^n* f^dr'r'^eyi^[-l3V^a{r)y ^' 

where n* = N^/r^ exp{~2f3Ec). 
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However, Eq. (Q) is more suitable for the derivation of the corresponding 
generahzed Kosterhtz equation. Employing the fugacity 



= (6) 

and the logarithmic length scale I = ln(r/ro) as defined in the BKT theory, we 
obtain, by taking the logarithm of Eq. (H) and differentiating with respect to I, 



dl 



Defining the dielectric constant as the ratio of bare and effective forces, 

_ dV/dr _ dV/dl 
^ ~ dV,s/dr ~ dV,s/dV 



(8) 



assuming a logarithmic form of the potential, V ^ ln(r/ro) = q^l, and defin- 
ing the stiffness constant of BKT theory by 

0a'^ 

we eventually find 

y' = (2-7rA>-27rV- (10) 

This equation contains the additional term —2Tr^y^ as compared with the orig- 
inal Kosterlitz equation. The corresponding initial condition, 

yo = yil = 0)=NoeM~PE,), (11) 

is the same as in BKT theory since the geometric correction results in a factor 
of unity in Eq. for the smallest pairs. 

The additional term in Eq. (|lo| ) is proportional to y^ as predicted by Halpe- 
rin 1^ . ft is also of the same form and has the same sign as the term found by 
Amit et al. ||lo| . Note that the notation of Ref. differs from the usual one of 
BKT theorypj. The implications of this term are discussed in Sec. 0. Here, we 
only note that since y > by definition and y' < for sufficiently large y, the 
fugacity y is bounded for any K, whereas y diverges for T > Tc in the original 
BKT theory. Equation ^ then implies that the density v falls off at least as 
fast as for large r so that the total density n also remains finite. 

3 Polarization and local field correction 

The question arises of whether higher order terms have also to be included in 
the second Kosterlitz equation in order to consistently describe higher vortex 
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densities. To address this question we first rederive the equation for the stiffness 
coefficient K starting from linear response theory and then discuss the necessary 
changes due to finite vortex densities. In the present section we use the language 
of the Coulomb gas model, i.e., we view the vortices as charged hard-core par- 
ticles interacting via the two-dimensional Coulomb potential V = ln(r/ro). 

In linear response theory the dielectric function e(k) = l/(k)/T4ff(k) is de- 
rived, 

' - ^^^(k), (12) 



E(k) 



= 1 - 



fc2 



where (?(k) is the Fourier transform of the charge (vorticity) density correlation 
function, 

5(r) = (p(O)p(r)). (13) 



The macroscopic dielectric constant is defined as 

— = lim — 
Coo /^^o e(k) 



(14) 



Assuming that l/foo is finite, we see from Eq. ( p2[ ) that limj.^o and lim^-^o dg(]<.)/dk 
have to vanish so that jl^, |l^ 



1 



k^O 



1 -7r/3 / d^rg{r) 



2 



/■OO /'27T 

l + TT/S drr d<p g{r)r'"' cos^ (j) 
Jo Jo 

/>oo 

1 + TT^/? / drr^g{r). 
Jo 



Here, we always have 1/eoo < 1. 

One now defines a scale-dependent dielectric constant ||l^, , 



(15) 



e(r) 



dr' r'^g{r'), 



(16) 



which includes the contribution of pairs of sizes r' < r only. Since we employ 
a hard-disk model we have introduced a lower cut-off tq in the integral. Note 
that e(r) is not the Fourier transform of e(k) as given by Eq. (p^). 

Since the correlation function g is not known exactly we need an approx- 
imation for this quantity. As discussed above, we can always decompose the 
vortex system into neutral pairs. We expect the correlation function to consist 
of three terms: (a) A ^-function term denoting the autocorrelation. This term 
is irrelevant because of the lower cut-off vq. (b) A term from the partner of 
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the test vortex, which should be proportional to the pair size distribution i^ir). 
(c) A term from the polarization of the other pairs. Here, we make the usual 
assumption that only smaller pairs are polarized in the field of a large one |^ . 
Neglecting the autocorrelation term, we write 

9{r) = 5pair(r) + 5poi(f)- (17) 
The relation between (/pair and the pair density v is given by 

9,Mr) - (18) 
7rr 

taking into account that v has already been integrated over the angle. To de- 
termine ^poi we employ a dielectric approximation. We assume that the smaller 
pairs form a dielectric gas characterized by the dielectric constant e(r), which is 
to be obtained self-consistently. The appearance of e(r) instead of too is a result 
of the approximation that only pairs smaller than r contribute to the screening. 

If two test charges, g at and —q at r', are brought into a neutral medium 
with dielectric constant e{r'), the total charge density, including the polarization 
charges, is given by 

P*W = ^^(r)-^<5(r-r'). (19) 

In the present approach the effect of screening is to reduce the charges of a pair 
of size r to the effective values ±g/e(r). Thus, the effect of screening on the 
correlation function can be expressed as 

g{r) = gpah.(r) + ffpoi(r) = -^-^TT' (^0) 

7rre^(r) 



From Eq. ( |l6| ) we get 



-L = l-TrPq' f^dr'r'^^. (21) 



Employing the definitions (^) and (|^), we obtain 

lF^'-|?i''"'''<''>^""'' '''' 

and, after differentiating with respect to Z, 

K' = -iTr^y^K^. (23) 

This equation is identical to the one obtained in BKT theory. 

In the above considerations the smaller vortex-antivortex pairs are treated as 
a continuous polarizable medium. In particular, local field effects are neglected 



7 



in the correlation function given by Eq. (pO|). Consequently, the dielectric con- 
stant of Eq. ( pi] ) does not include these effects, either. Here, we employ a 
Clausius-Mossotti type formula to obtain a better approximation for e. In two 
dimensions the relation between the dielectric constant neglecting local field 
effects, Ea, and the full dielectric constant e is given by 



(24) 



Taking from Eq. (|6|) we thus make the ansatz 

e{r) ^ 1 + f /5 J^^dr' r'^g{r') 

l-^l3q^ J^^dr' r'^u{r')/e^{r') 
1 - ^ Pq"^ dr' r''^iy{r')/e^{r') 



(25) 



where g[r) is not the (unknown) exact correlation function, but the one with- 
out local field corrections, Eq. (20). Inserting the definitions (||) and (||) and 
expanding for small y, we obtain 



K 



27T 



di'y'inK'in 



dl' y\l')K\l') 



Differentiation yields 

K' = -Aw^y^K^ 



1 + 



8^ 



dl'y'(l')K^l') 



(26) 



(27) 



This expression is similar to the predicted one [|| ^: The second term on the 
right hand side is of fourth order in y. However, it is not simply proportional to 
y'^. Note that this result is qualitatively different from the one of Amit et al. pO| . 
These authors effectively perform an expansion for a different small parameter, 
namely K — 2/tt instead of y, which makes direct comparison difficult. 

To second order in y, the term in brackets in Eq. ( p7|) is given by 2~2ttK / [3q^ ^ 
as can be seen from Eq. (^). Thus we find, to fourth order in ?/, 



K' = -An^y^K^ 



2ttK 



(3q^ ) 



(28) 



For pairs of minimum size rp there is no screening (e = 1) and, therefore, the 
initial condition reads 

Ko^K{l = 0) = ^. (29) 
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Equation ( p8[ ) can thus be rewritten as 

K' = -An^y^K^ ^2 - . (30) 

4 Discussion of the generalized equations 

The generahzed Kosterhtz equations now have the form 

K' = -ATT^y^K^(^2^^y (31) 

y' = (2-7ri^)y-27rV, (32) 

together with the initial conditions given by Eqs. (|ll|) and (H). We proceed 
to discuss these equations. Note that the equation for K' contains the initial 
value Kq = I2itT as an additional parameter so that the derivatives at a point 
(X, y) are not exclusively determined by K and y, as opposed to the original 
BKT theory. 

We are interested in the behaviour of Eqs. (|3l|) and ( p2| ) on large length 
scales, I oo. The equations have an attractive fixed point at = 0, y = I/tt 
and an attractive line of fixed points at K > 2/tt, y = 0. The first fixed point 
corresponds to a screened interaction at large distances {K — > 0); it is the high- 
temperature fixed point. Contrary to BKT theory, however, where this fixed 
point is a,t K = 0, y = oo, the fugacity remains finite. We note here that our 
approach, which essentially treats all vortices as bound in pairs and thus has 
no free vortices, predicts the correct, metallic form of screening for the high- 
temperature phase. Of course one has to check whether the decomposition of 
the vortex system into pairs is justified for T > Tc- This is done below. The 
line of fixed points corresponds to a vanishing density of large pairs {y 0) 
together with a logarithmic interaction at large distances {K remaining finite). 
This is the bound-pair phase. 

In BKT theory analytical expressions for the trajectories, i.e. for the curves 



{K{l),y{l)), can be found. The same is not true for Eqs. (31 ) and (|32|). However, 
these equations are not difficult to integrate numerically. Since the equations 
contain the additional parameter Kq, every initial point now defines its own 
trajectory. In Fig. |l| trajectories for Ec/q^ — 0.5 and several temperatures (thin 
solid lines) as well as for Ec/q^ = 0.2 and several temperatures (dashed lines) 
are plotted. The dotted lines denote curves of initial values for Ec/q^ — 0.2 and 
Ec/q^ = 0.5, respectively, and varying T jq^ . The trajectories are shown in K-y 
space instead of the more usual x-y space, where x = — 1 + 2/'kK, to empha- 
size the high-temperature behaviour. The trajectories in x-y representation are 
similar to the ones of standard BKT theory. 

An important question is which initial values belong to which phase, or put 
mathematically, which points belong to the basin of attraction of the high-tem- 
perature fixed point and of the line of low-temperature fixed points, respectively. 
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As in the original BKT theory, the two basins of attraction are separated by 
the initial values which flow to X 2/7r, ?/ = in the limit I oo. These 
points form the so-called separatrix, which can be obtained numerically and is 
shown in Fig. |l| as the thick solid line. The separatrix is found to lie below 
the BKT result, i.e., the high-temperature phase is stabilized. This result is 
not surprising since local field effects tend to increase e, thereby weakening the 
renormalized vortex interaction. 

Although the trajectories themselves cannot be obtained analytically, limit- 
ing forms of the separatrix and the leading temperature dependence of Koo = 
K{1 oo) for Kq 2/7r and y ^ 1 can be deduced from the recursion relations 
and (^^. First, we consider the renormalized stiffness coefhcient i^oo- For 
T > Tc the flow is to the high-temperature flxed point and we have Koo — 0. 
At the transition we recover the universal value Koc = 2/7r of BKT theory. To 
find the behaviour for T < Tc, we consider Eqs. ( ^) and (|2|) to leading order 
only. 



K' = -Wny^, 



iy 



2\l 



= -2ny^ K 



(33) 
(34) 



These equations lead to 



K'- 



'2ttK' K 



(35) 



with the initial conditions K{{)) = Pq^ /2tt and K'{Q) = The solution 

is 

o /7-<" / r rzr / o \ i \ 

(36) 




with C = 7r{Ka — 2/tt) — IGiry^. In the limit Z — > oo we thus obtain 




Ko-- -Wy'o- 

TT 



(37) 



We know from the analysis of the fixed points that Koo = 2/tt a±T — Tc- There- 
fore, to the present order the radicand in Eq. (^) vanishes at the transition. 
To the same order we thus find 



4 V TT 



(38) 



as the limiting form of the separatrix, which is the same result as given by BKT 
theory. 
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From the initial conditions (|Tl|) and ( p9| ) we see that the leading temperature 
dependence of the radicand in Eq. (B^) is linear so that 




^oo--+4Ws^^ (39) 



for small Tc — T > 0. Here, B is some function of the core energy independent 
of temperature. Equation ( ^9|) also has the same form as in BKT theory 0|. 
The typical square-root cusp and the universal jump are thus left unchanged by 
the higher order terms. 

The same is true for the typical length scale which is interpreted as the 
correlation length of the superconducting order parameter. This length is also 
used as the aforementioned cut-off in BKT theory. Here, it does not have such 
a significance and we only discuss it for completeness. We define this quantity 
as the length scale on which the trajectories start to be drawn towards the high- 
temperature fixed point, i.e., we define K{r = ) = 2/7r. The length scale ^-|- 
thus diverges for T ^ Tc from above. For T < Tc the high-temperature fixed 
point never comes into play and we have ^+ = 00. Under the same assumptions 
as above and starting again from Eqs. ( ^3| ) and ( ]3^ ) we find 

In ^ ^ ^ . (40) 

ro ^^16%' - (Ko - 2/^)2 

Since the radicand is again linear in T to leading order, we reobtain the BKT 
result 

4, 



«*"^"°"=l7f^j (41) 

for small T — Tc > with b independent of temperature. 

The stability of the well-established predictions for Koc and ^+ in the pres- 
ence of higher order terms lends additional support to the general concept of 
BKT. As far as these quantities are concerned, however, both the original the- 
ory and the extension presented here are in agreement with experiment. Thus, 
experiments such as measurements of the current-voltage characteristics of su- 
perconducting films, which effectively measure K, cannot distinguish between 
the two approaches. 

Such a distinction may be possible, however, using the total pair density n. 
In Sec. ^ we already mentioned that y is bounded in the present case and that, 
therefore, n must remain finite. We now discuss this point in more detail. 

From Eqs. (|l|) and (^) the total density in natural units is given by 

/•oo 

nrl =2tt dle~^^y^{l). (42) 
^0 

This integral cannot be evaluated analytically since y^(Z) is unknown. One can. 
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however, obtain a rigorous statement concerning the continuity of n. Substitu- 
tion of u = e^' for / yields 

nrl ^2tt [ duuy'^{u). (43) 



^0 

In this expression, is a continuous function of I, and thus of u. Furthermore, 
is a smooth, i.e. arbitrarily often differentiable, function of the initial values 
Kq and yo for any finite /, see Fig.0. Thus y'^ is a smooth function of temperature 
for finite I (u > 0). On the other hand, y^ is not a smooth, or even continuous, 
function of T for I — oo (u = 0) since it has a jump at Tc- The integrand in 
Eq. (|4^), however, is smooth in T even for u = 0. Thus we have an integral over 
a finite interval over a function that is continuous in the integration variable 
and smooth in the external parameter T. Under these suppositions the integral 
nrp is a smooth function of temperature even at the transition. This behaviour 
is connected with the observation that the phase transition is of infinite order 

i- . 

This result is in agreement with BKT theory, where, however, one has to 
introduce a cut-off ^+ to obtain it. Since the definition of ^+ is somewhat 
arbitrary, as discussed in Sec. |^, the argument given here appears to be based 
on firmer grounds. 

To obtain numerical results for the total density n, the generalized Kosterlitz 
equations (^) and (|2h are integrated numerically up to Z = 20, corresponding 
to a system size of er^ro, for various values of T/q^ and Ec/q^- We assume 
A*o = 1. The results for y are used to approximate the integral in Eq. (|4|). The 
results are plotted in Fig. ^. As expected, the density increases with temperature 
and decreases with increasing core energy. Note that the approach presented 
here is not limited to T/q^ < 1/4, where 1/4 is the maximum value of Tc/q'^. 
The density starts to increase below Tc, but the larger part of the increase 
takes place above Tc- As expected, there is no feature at Tc (open circles in 
Fig. H). For high temperatures the density saturates. This saturation and 
in fact the convergence of the density for T > Tc are due to the geometric 
correction of Sec. The temperature dependence of the total density is, in 
principle, experimentally accessible through, e.g., flux noise and NQR 
measurements on superconducting films. 

One may ask whether the theory is still applicable for T > Tc, where the 
average pair size approaches the typical distance between neighboring pairs. If 
both quantities are of the same order the decomposition of the vortex system 
into small pairs becomes meaningless. To address this issue we compare the 
pair separation rp with the average pair size (r). The pair separation is just 
Tp — l/^/n. The pair size is 

1 rdle-'y^m 
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This quantity can be calculated similarly to the pair density. The result for 
Ec/q^ = 0.5 is depicted in Fig. |3[ The important point is that (r) (dashed line) 
remains small compared to Tp (solid line) even in a significant temperature range 
above Tc- In this regime large pairs are "broken" in the sense that the interaction 
takes on a screened form for large distances, but nevertheless the pairs do not 
overlap considerably and the decomposition of the vortex system into pairs is 
justified. Therefore, the theory is applicable in the proximity of T^, in addition 
to the low-temperature phase. The inset of Fig. || shows the corresponding 
quantities for E^jq^ = 0.05. One sees that the range of applicability for T > Tc 
is narrower in this case, but it still exists. 

We thus have reason to believe that the extended theory is valid in a temper- 
ature interval around Tc for arbitrary core energy Ec, and, therefore, that the 
transition is of infinite order for any Ec- This result contradicts Refs. [|l3[ [l5|, 
which predict a first-order transition for small Ec/q^- However, we cannot ex- 
clude the possibility that terms of even higher order in the Kosterlitz equations 
bring about such a transition or that the expansion for small y does not work 
at all for small Ec/q^- 

To conclude, straightforward physical models have been presented for higher 
order terms in the Kosterlitz equations. These terms are of a form similar, but 
not identical, to earlier predictions ||, || . A geometric correction, which ensures 
the correct enumeration of pairs, leads to the disappearance of the unphysical 
divergence of the total pair density for T > Tc- Starting from an exact ex- 
pression for the dielectric constant, the polarization of the vortex system and 
local field effects are taken into account, leading to a stabilization of the high- 
temperature phase. The behaviour of the stiffness constant and the correlation 
length close to Tc is identical to the predictions of BKT theory so that ex- 
periments on these quantities, e.g. measurements of the nonlinear resistance of 
superconducting films, cannot distinguish between the two approaches. The 
employment of the Clausius-Mossotti formula in Sec. § does not lead to easily 
testable predictions. However, it is comforting to know that the next higher 
order term in the equation for K does not destroy the BKT behaviour. The 
enumeration rule of Sec. |^, on the other hand, leads to new quantitative predic- 
tions for the total pair density even above Tc, which can, in principle, be tested 
by flux noise and NQR techniques. In short, it has been shown that the BKT 
approach can be extended in such a way that the temperature regimes below 
and above Tc are treated on equal footing. 
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Figures 



Figure 1: Trajectories of the generalized Kosterlitz equations ( |3l| ) and ( p2[ ) 
for Ec/q^ = 0.5 and varying T /q^ (thin soUd Hues), and for Ec/^ — 0.2 and 
varying T/g^ (dashed Hues). The separatrix is denoted by the thick sohd Hne, 
the bound-pair phase is below that line. The dotted lines denote initial values 
{I = 0) for Ec/q^ — 0.2 and 0.5, respectively, and varying temperature. 



Figure 2: The total pair density n as a function of T/q^ for Ec/q^ ~ 0.02, 0.05, 
0.1, 0.2, 0.5, 1 (from left to right). The open circles denote the respective values 
of Tc/q^ and the corresponding densities. 



Figure 3: The average pair size (r) (dashed line) and the average separation be- 
tween neighboring pairs (solid line) for Ec/q^ = 0.5. The filled triangle denotes 
Tc- The inset shows the same quantities for Ec/q^ = 0.05. 
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